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ABSTRACT 

Apollonian circle packings arise by repeatedly filling the interstices between four mutually 
tangent circles with further tangent circles. Such packings can be described in terms of the 
Descartes configurations they contain, where a Descartes configuration is a set of four mutually 
tangent circles in the Riemann sphere, having disjoint interiors. Part I showed there exists 
a discrete group, the Apollonian group, acting on a parameter space of (ordered, oriented) 
Descartes configurations, such that the Descartes configurations in a packing formed an orbit 
under the action of this group. It is observed there exist infinitely many types of integral 
Apollonian packings in which all circles had integer curvatures, with the integral structure 
being related to the integral nature of the Apollonian group. Here we consider the action of 
a larger discrete group, the super-Apollonian group, also having an integral structure, whose 
orbits describe the Descartes quadruples of a geometric object we call a super-packing. The 
circles in a super-packing never cross each other but are nested to an arbitrary depth. Certain 
Apollonian packings and super-packings are strongly integral in the sense that the curvatures 
of all circles are integral and the curvature x centers of all circles are integral. We show that (up 
to scale) there are exactly 8 different (geometric) strongly integral super-packings, and that 
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each contains a copy of every integral Apollonian circle packing (also up to scale). We show 
that the super-Apollonian group has finite volume in the group of all automorphisms of the 
parameter space of Descartes configurations, which is isomorphic to the Lorentz group 0(3, 1). 

Keywords: Circle packings, Apollonian circles, Diophantine equations, Lorentz group, Coxeter 
group 

1. Introduction 

Apollonian circle packings are arrangements of tangent circles that arise by repeatedly filling 
the interstices between four mutually tangent circles with further tangent circles. A set of four 
mutually tangent circles is called a Descartes configuration. Part I studied Apollonian circle 
packings in terms of the set of Descartes configurations that they contain. It is observed that 
there exist Apollonian circle packings that have a very strong integral structure, with all circles 
in the packing having integer curvatures, and rational centers, such that curvature x center is 
an integer vector. We termed these strongly integral Apollonian circle packings. An example 
is the (0,0, 1, 1) packing pictured in Figure ^ with the two circles of radius 1 touching at the 
origin, and with two straight lines parallel to the x-axis. 




Figure 1: The integer Apollonian packing (0,0, 1, 1) 

Part I gave an explanation for the existence of such integral structures. This uses a co- 
ordinate system for describing all (ordered, oriented) Descartes configurations T> in terms of 
their curvatures and centers, which forms a 4 x 3 curvature- center coordinate matrix Mp, 
and a more detailed coordinate system, augmented curvature- center coordinates, using 4x4 
matrices W-p. The strongly integral property of a single Descartes configuration is encoded 
in the integrality of the matrix . The set of all (geometric) Descartes configurations in an 
Apollonian packing can be described as a single orbit of a certain discrete group A of 4 x 4 
integer matrices; algebraically there are 48 orbits of ordered, oriented Descartes configurations 
giving rise to the same geometric packing, which correspond to the 48 possible ways of ordering 
the circles and totally orienting the configuration. The integrality of the members of A is the 
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source of the strong integrality of some Apollonian circle packings. As a consequence of this 
group action, if a single Descartes configuration in the packing is strongly integral, then they 
all are; hence every individual circle in the packing is strongly integral. 

There are infinitely many distinct integral Apollonian circle packings. Two more of them 
are pictured in Figures 2 and 3, the (— 1, 2, 2, 3)-packing and the (— 6, 11, 14, 15)-packing, re- 
spectively. Any integral packing can be moved by a Euclidean motion so as to be strongly 
integral, as will follow from results in this paper. 




Figure 2: The integer Apollonian packing (—1,2,2,3) 



Part I introduced some further integral actions on Descartes configurations, involving a 
"duality" operator D leading to a dual Apollonian group A^~. Combining this group with the 
Apollonian group leads to a large group of integer 4x4 matrices, the super- Apollonian group 
A s , which acts on the set of all (ordered, oriented) Descartes configurations. Part I showed 
that the super- Apollonian group is a hyperbolic Coxeter group. It defined an Apollonian super- 
packing to be an orbit .A 5 [2?] of a single Descartes configuration T> under this group. Such a 
super-packing is called integral if the initial Descartes configuration has integer curvatures, and 
is strongly integral if the curvature-center coordinates Mp of the initial Descartes configuration 
are integral. These properties hold for all Descartes configurations in the packing if they hold 
for one. 

In this paper we study the geometric structure of super-packings, and their integrality 
properties. A striking geometric fact is that the circles in a super-packing do not overlap, as 
shown in §3. Figure 0] shows circles of curvature at most 65 in the super-packing generated 
from the (0,0,1,1) configuration in Figure Q above. (The generating Descartes configuration 
is indicated with slightly darker lines.) Here the horizontal and vertical scales of the figure are 
roughly from —2.2 < x < 2.2 and —2.2 < y < 2.2. This picture is representative of the whole 
super-packing; we show in §6 that this super-packing is periodic under the two-dimensional 
lattice generated by the shifts x — > x + 2 and y — > y + 2. 
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We cannot picture the whole super-packing because the circles in it are dense in the plane; 
they nest inside each other to an arbitrary depth. One can show that the set of points lying in 
the interior of an infinite nested sequence of such circles has full Lebesgue measure. However 
there is interesting structure visible if we turn up the scale of the circles. In §6 in Figure |H1 we 
show all circles of curvature at most 200 inside the super-packing above, inside the unit square 
< x < 1, <y < 1. 

This paper is concerned with integrality properties of some super-packings, and we study 
successively stronger integrality properties. In §4 we require the curvatures to be integral, in 
§5 and §6 both the curvatures and curvature x center must be integral, and in §8 the augmented 
curvature-center coordinates of all Descartes configurations in the super-packing must be inte- 
gral. At this last level we show that there are exactly 14 such super-packings, viewed as rigid 
geometric objects. In §2 we give a more detailed description of our results. 

The general framework of this paper was developed by the second author (JCL), who also 
did much of the writing. This paper is an extensive revision of a preprint written in 2000, and 
adds new results in §3 and §6. 

Acknowledgments. The authors are grateful for helpful comments from Andrew Odlyzko, 
Eric Rains, Jim Reeds and Neil Sloane during this work. They also thank the reviewer for 
incisive comments. 
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Figure 4: An initial part of the (0,0, 1, 1) super-packing (square of sidelength 4.4) 
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2. Summary of Main Results 



We consider Apollonian super-packings. Analogously to part I, an Apollonian super-packing 
may be considered as either an geometric object or an algebraic object, as follows. 

(i) [Geometric] A geometric Apollonian super-packing is a point set on the Riemann sphere 
C = M? U {00} , consisting of all the circles in four orbits of a certain group Gj^s{T>) of Mobius 
transformations inside the conformal group M6b(2) acting on the four circles {C\, C2, C3, C4} 
in a given Descartes configuration T>. The group G^s(T>) depends on T>. 

(ii) [Algebraic] An ( algebraic ) Apollonian super-packing is a set of ordered, oriented Descartes 
configurations, given by 48 orbits of the super- Apollonian group .A 5 [2?]. The augmented 
curvature-center coordinates of its elements are A s Wx> ■= {UW-p : U G A s }. 

A geometric super-packing can be described in terms of its unordered, unoriented Descartes 
configurations. From this viewpoint, each geometric Apollonian super-packing corresponds to 
48 different algebraic super-packings; there are 24 choices of ordering the four circles and 2 
choices of total orientation of the configuration. We can consider geometric super-packings as 
unions of a countable number of Apollonian packings. This leads to interesting questions con- 
cerning the way these Apollonian packings are embedded inside the geometric super-packing. 
We note that as a point set, the geometric super-packing is invariant under the group action 
G^s (T>). However it is not a closed set, and its closure is the entire Riemann sphere R 2 U {00}. 

A large part of the paper considers integrality properties of curvatures and centers of some 
super-packings. These questions are mainly studied using algebraic super-packings, although 
we also consider questions concerning the associated geometric super-packing, such as its group 
of symmetries under Euclidean motions. 

In §3 Theorem 13.11 shows that each geometric super-packing is a packing in the sense that 
the circles in it do not cross each other transversally, as mentioned above. Theorem l3.2l specifies 
certain subcollections of geometric super-packings which are genuine packings in the sense that 
the interiors of the circles do not overlap. 

In §4 we study integer super-packings, in which all circles have integer curvatures. Integer 
super-packings are classified up to Euclidean motions by a single invariant, their divisor g, 
which is the greatest common divisor of the curvatures in any Descartes configuration in the 
packing. Theorem 14 . 1 1 shows that for each g > 1 there is a unique such integral super-packing, 
up to a Euclidean motion. We also show in Theorem 14.31 that for each geometric Apollonian 
circle packing that is integral, there exists a Euclidean motion taking it to one that is strongly 
integral. 

In §5 we study strongly integral super-packings, which are those whose curvatures are 
integral and whose curvature x center is also integral. Strongly- integral super-packings are 
geometrically rigid: Theorem 15.11 shows that for each integer g > 1 there are exactly 8 strongly 
integral geometric super-packings which have divisor g. Here we do not allow the packings to 
be moved by Euclidean motions. 

In §6 we study the relations between integer Apollonian packings and strongly integral 
super-packings. Without loss of generality we restrict to primitive integer super-packings, 
those with divisor 1. Theorem 16 . II shows that each of the 8 kinds of these has a large group of 
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internal symmetries, forming a crystallographic group. For convenience we fix one of them and 
call it the standard strongly integral super-packing; results proved for it have analogues for the 
other seven. Theorem 16.21 shows that each primitive integral Descartes configuration (except 
for the (0,0,1,1) configuration) occurs in this packing with the center of its largest circle being 
contained in the closed unit square < x < 1, < y < 1, and the location of this circle center 
is unique. Theorem 16.31 deduces that the geometric standard strongly integral super-packing 
contains a unique copy of each primitive integral Apollonian packing, except for the (0, 0, 1, 1) 
packing, having the property that the center of its bounding circle lies inside the closed unit 
square. Figures El 101 an d picture the locations of all primitive integer Apollonian packing 
with bounding circle of curvatures 6, 8 and 9, respectively. The unit square is indicated by 
slightly darker shading in the figures. Note that in Figure El the three Apollonian packings are 
generated by Descartes configurations with curvature vectors (—6,7,42,43), (—6,10,15,19) 
and (—6, 11, 14, 15); these are root quadruples in the sense of Sec. 4]. 




Figure 5: Integer Apollonian packings with bounding circle of curvature 6. 

In §7 Theorem 17. ll shows that the super- Apollonian group is a finite index normal subgroup 
of the group Aut{Qo^) of integral automorphs of the Descartes quadratic form. The latter 
group can be identified with an index 2 subgroup of the integer Lorentz group 0(3, 1,Z), and 
this identification allows us to identify the super-Apollonian group with a particular normal 
subgroup A s of index 96 in 0(3, 1,Z), defined after Theorem 17. II 

In §8 we study super-packings all of whose Descartes configurations are super-integral in 
the sense that their augmented curvature-center coordinate matrices are integer matrices. 
Theorem 18. II shows there are exactly 14 geometric super-packings of this kind. 

The last section §9 makes a few concluding remarks. 

The Appendix considers minimal conditions to guarantee that a Descartes configuration is 
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Figure 6: Integer Apollonian packings with bounding circle of curvature 8. 




Figure 7: Integer Apollonian packings with bounding circle of curvature 9. 



7 



strongly integral. Theorem 1 1 U . 1 1 shows that a configuration is strongly integral if (and only if) 
three of its four circles are strongly integral. 

3. Geometric Apollonian Super-Packings 

In this section we consider properties of geometric Apollonian circle packings. We view such 
a super-packing as a point set on the Riemann sphere C = C U {oo}. We first note that it 
is not a closed set. It is not hard to show that its closure is the whole Riemann sphere, 
Each geometric Apollonian packing has a group invariance property under a certain group of 
Mobius transformations which depends on the super-packing, which is the group generated by 
the countable set of groups of Mobius transformations that leave invariant some Apollonian 
packing contained in the super-packing. 

Our object is prove the following "packing" property of geometric super-packings. 

Theorem 3.1. A geometric Apollonian super-packing is a circle packing in the weak sense 
that no two circles belonging to it cross each other transver sally. Circles in the geometric 
super-packing may be nested, or tangent to each other. 

Before giving the proof, we describe the nature of the geometric packing in terms of 
nesting of circles. We view the packing ,4 s [2?q] as generated from an initial (positively ori- 
ented) Descartes configuration Vq, by multiplication by a finite set of generators of the super- 
Apollonian group. Each circle in the super-packing has a well-defined nesting depth d (relative 
to the generating configuration Dq) which counts the number of circles in the packing which 
include C in their interior. Here the notion of "interior" is defined with respect to the initial 
Descartes configuration T>q. The Apollonian group generators move "horizontally", leaving 
constant the nesting depth of any circles they produce. The dual Apollonian group generators 
move "vertically", by reflecting three of the circles in a configuration into the interior of the 
fourth circle, they increase the nesting depth by one. We show there is a unique "normal form" 
word of minimal length in the generators that produces a Descartes configuration T> containing 
C. The nesting depth of C is exactly equal to the number d of generators of A ± that appear 
in this normal form word. The circles at nesting depth are those circles in the Apollonian 
packing generated by T>q. Each circle C at nesting depth k > 1 contains a unique Apollonian 
packing, consisting of it plus all circles at depth k + 1 nested inside it. 

Proof of Theorem 13.11 We view the geometric super-packing on the Riemann sphere C = 
C U {oo}, so the initial Descartes configuration T>q consists of four circles on the sphere. In 
this case each circle defines a spherical cap. We choose an ordering and orientation of T>o (this 
does not affect the geometry), requiring that T>q have positive (total) orientation. 

Let the super-packing be A [T>o] = A s Wx> . We consider the effect of the super- Apollonian 
group generators acting on the left on the matrix Wp . The Descartes configurations in the 
super-packing are given by V = U m U m _i • • • Ui[2?o] with each S {Si, S2, S3, S4, S^, S%, S3-, S4-}. 
We consider words in the group measured by their length m. The stage m circles will consist 
of all new circles added using products of m generators. We may without loss of generality 
restrict to normal form words, which are those that satisfy the two conditions: 

(i) If U fc = Sj, then U fc+ i ^ Si and U fc+ i ^ Sf with j / i. 
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(iiJIfU^S-L.thenUfc+i^S^. 

Equivalently, looking backwards, if Uk+i = then either = Sj or else = Sj- for 
some j 7^ i. A word may be put in normal form by canceling adjacent equal generators, since 
all S 2 = (S^) 2 = I, and by moving towards the right 2 in the word as far as possible any 
generator S^~, using the property that it commutes with all Sj with j / i. These operations 
eventually put a word in normal form, without increasing its length. The operations do not 
change the Descartes configuration T> it represents. 

We prove the theorem by induction on the number of symbols m in a normal form word, 
which we call the stage of the induction. The induction hypotheses at stage m are as follows. 
Here we let Cj refer to the circle at row i of the associated (ordered, oriented) Descartes 
configuration. 

(1) Each normal form word of length m produces either one or three new circles, according 
as U m = Sj, where it is the circle Cj or U m = Sj-, where it is the three circles Cj with i 7^ j. 

(2) The nesting depth of any new circle produced at this stage is equal to the number of 
occurrences of a letter in A 1 " in its generating normal form word. 

(3) Each new circle produced has empty interior, when it first appears. 

In particular hypothesis (3) implies that all circles produced at stage m have disjoint 
interiors, and that each such circle contains no circles from earlier stages in its interior. If 
the induction is proved, then hypothesis (3) guarantees that the nesting depth of a circle is 
well-defined when it is first produced, because no new circle will ever include it in its interior. 
Hypothesis (3) also guarantees that all circles produced by the end of stage m + 1 do not cross. 
As a consequence no two circles in the packing cross, using property (3) applied at that level 
m which is the greater of the levels of the two circles. Thus the theorem will follow. 

The base case m = of the induction is immediate, consisting of the initial Descartes 
configuration Vq. We now show the induction step for m + 1, given m. We are given a 
Descartes configuration V' = U m +i[P] = U m+ iU m • • • Ui[2?o] with a normal form word. 

To establish hypothesis (1) for m + 1, suppose first that U m +i = Sj~. We assert that the 
2-th circle Cj of T> was a new circle produced at stage m. For either U m = Si, in which 
case it was the unique new circle in V by induction hypothesis (1) at stage m, or U m = Sj~ 
with j / i, in which case it was one of three new circles produced at stage m. By induction 
hypothesis (2) the circle Cj has empty interior at stage m. The three circles {Cj : j 7^ i} in 
the new configuration D' are contained in the interior of Cj are therefore new circles. They do 
not cross, being part of a Descartes configuration. Thus hypothesis (1) holds for m + 1 in this 
case. 

Suppose next that U m+ i = Sj, so the possible new circle is C i . If T>k := UfcUfc_i • • • Ui[X>o] 
is the maximal length subword such that = Sj- for some j, with 1 < k < m then V belongs 
to the Apollonian packing generated by T>k, since all subsequent generators belong to the 
Apollonian group. If no such k exists, then V is in the Apollonian packing generated by the 
original Descartes configuration T>q. For k > 1, this Apollonian packing is entirely contained 
in the interior of a bounding circle C = Cj first produced at stage k — 1. At that time the 
interior of C was empty, by induction hypotheses (2) and (3) at stage k— 1. The only Descartes 
configurations that ever can enter the interior of the circle C, must do so by a reflection in C 
and these are exactly those normal form words starting with initial segment UfcU/%_i • • -Ui. 

2 That is, moving it towards the beginning of the word. 
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(This follows from uniqueness of a circle when it is created, hypothesis (3) applied at stage k.) 
The words contain this initial segment at the same depth, with all subsequent letters in A, fill 
out an Apollonian packing at this depth. In particular each such normal form word produces 
one new circle in this Apollonian packing. Recall that all the circles in the Apollonian packing 
inside the bounding circle C have disjoint interiors (Theorem 4.1 of part I). These circles 
all have the same depth, and (1) holds in this case. Normal form words that have another 
subsequent generator in A 1 - confine the resulting Descartes configuration to the inside of a 
single circle in the Apollonian packing already produced, and all longer words with this 

prefix are entirely contained inside this circle. In particular, they do not coincide with C i . the 
new circle produced by the normal form word corresponding to V. It follows that the circle 
C i is new, so hypothesis (1) holds for m + 1 in this case. 

There remains the case where U m +i = Sj and all previous £ A. Then U m +i • • • Ui [Dq] 
belongs to the Apollonian packing generated by T>q, and is at depth 0. Here C i is new since 
the generation of an Apollonian packing creates one new circle at each step, and any word 
that contains an element of A 1 " moves the Descartes configuration inside an older circle in this 
packing, from which it cannot escape. Thus ^ is a new circle in this case, and hypothesis (1) 
holds. 

Hypothesis (2) holds for m + 1 in the case U m = Sj- because the three new circles produced 
have nesting depth one greater than the nesting depth of C\ at the previous level, to which 
induction hypothesis (2) applied. It also holds in the remaining case U m = Sj because the 
argument above showed that the nesting depth did not increase in this case. 

Hypothesis (3) holds if U m+ i = S^~. Now circle Ci was first created only at stage m, and 
the only other possible sequences leading to a Descartes configuration including this circle must 
start from V m and use a generator U m+ i = Sj with j ^ i. In all other cases the resulting 
Descartes configuration includes no circle inside Ci, so the interiors of the three new circles 
produced are empty at the end of stage m + 1. 

In the remaining case U m +i = Sj, we have already observed that the new circle Cj pro- 
duced, is disjoint from all other circles in the Apollonian packing -4[X\-] created by level m, 
which are necessarily contained in the bounding circle C of T>^. As mentioned above, all 
Descartes configurations containing a circle inside C must have an initial segment of their 
generating word, giving the unique normal form word that first generates C, at stage k < m. 
Other depth m + 1 words with this initial segment, and with all subsequent Uj drawn from 
the Apollonian group generators produce new circles in the Apollonian packing disjoint 
from Cj. Any normal form word at depth m + 1 with this initial segment, which contain some 
generator Sj~ after-wards, produce a Descartes configuration contained inside a circle of the 
Apollonian packing .AfX^] different from C . It follows that the interior of C[ is empty at the 
end of stage m + 1, as required. 

This completes the induction step and the proof. ■ 

Super-packings have some properties that are genuine packing properties. The proof of 
Theorem |23 established in hypothesis (3) shows that the finite set of circles at stage m of the 
construction, starting from a generating Descartes configuration D, all had disjoint interiors. 
It also gives the following stronger result. 

Theorem 3.2. For a geometric Apollonian super-packing given with a generating Descartes 
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configuration D, and each k > 1, i/ie sei o/ aZZ circles having nesting depth exactly k with 
respect to T> have pairwise disjoint interiors. These circles can be viewed as forming an infinite 
collection of Apollonian packings, each missing one circle; the missing circle is a bounding circle 
at depth k — 1. 

Proof. The proof of Theorem 13.11 shows that the nesting depth of a circle is well-defined. 
Circles at nesting depth k have disjoint interiors, since no two of them cross, and the only way 
for two of them to have an interior point in common is for one to be nested inside the other, 
which would violate the nesting ordering. 

Given a normal form word U := U m U m _i • • • Ui with U m = S^~, and containing exactly 
k elements of A 1 - = (S±, , S3, ), the set of all normal form words having V as prefix 
and all other letters JJj £ A for j > m produce all the circles in an Apollonian packing 
„4[D m ], with V m = XJ[T>o]. All these circles are at depth k except for the outer circle of V m , 
which is its i-th circle. Enumerating all possible such U as prefixes represents the set of nesting 
depth k circles as a collection of Apollonian packings, each excluding one circle when k > 1. ■ 

Remarks. (1) The proof of Theorem 13. II is a geometric analogue of the presentation for the 
super- Apollonian group proved in part I, [31 Theorem 6.1]. 

(2) The analogous result to Theorem l3 1 1 1 fails to hold in all dimensions n > 4, as explained in 
part III. The "nesting" property of the dual Apollonian configurations resulting from "vertical" 
moves still exists and works in all dimensions. However the "horizontal" motions moving 
spheres around in Apollonian packings produces spheres that cross in dimensions n > 4, see 

Lemma 4.1]. 

(3) We cannot easily visualize a geometric super-packing as a completed object because 
the circles in it are dense in the plane. We can however picture a partial version of it that 
pictures all circles of size above a given threshold, in some finite region of the packing. The 
integral super-packings we are most interested in have a periodic lattice of symmetries (see 
Theorem 16. ip . so it suffices to examine a finite region of the packing. Figure in §1 and 
Figure |H1 in §6 exhibit part of a super-packing. 

(4) Every circle C in a geometric super-packing _4 5 [D] has associated to it a unique Apollo- 
nian packing of which it is the bounding circle. If it is a depth k circle (relative to the starting 
configuration T>), then this Apollonian packing consists of all depth k + 1 circles contained in 
the interior of C. 

4. Integral Super-Packings 

An Apollonian super-packing is integral if it contains one (and hence all) Descartes configura- 
tion whose circles have integer curvatures. 

An invariant of an integral super-packing is its divisor g, which is the greatest common 
divisor of the curvatures of the circles in any Descartes configuration in the super-packing. 
The quantity g is well-defined independent of the Descartes configuration chosen, using the 
relation Mp/ = UMp between two such configurations, where U G A s . Since U is an integer 
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matrix with determinant ±1, it preserves the greatest common divisor of each column of the 
integer matrix Mp. Here the first column encodes the curvatures of the circle. 



Theorem 4.1. For each integer g > 1 there exists an integral Apollonian super-packing with 
divisor g. The associated geometric super-packing is unique, up to a Euclidean motion. 

As an immediate corollary of this result, the integral super-packing with divisor g contains 
at least one copy of every integral Descartes configuration having divisor g. For each such 
a Descartes configuration generates an integral super-packing with divisor g, so the corollary 
follows by the uniqueness assertion. 

We defer the proof of Theorem 14.11 to the end of the section. It is based on a reduction 
theory which finds inside any such integral super-packing a Descartes configuration having 
particularly simple curvatures. 

Theorem 4.2. LetD be an integral Descartes configuration, with divisor g := gcd(&i, 62, 63, 64). 
Then the Apollonian super-packing A S [T>] generated by T> contains a Descartes configuration 
having curvature vector a permutation of either (0,0, g,g) or (0, 0, — g, — g), with the former 
case occurring i/ &i + 62 + 63 + 64 > and the latter case if b\ + 62 + 63 + 64 < 0. 

Proof. Since the super-Apollonian group A s preserves the (total) orientation of Descartes 
configurations, it is sufficient to show that for a positively oriented integral Descartes config- 
uration T> with curvatures (bi, b2, 63, 64), b\ + 62 + &3 + &4 > 0, there exists U £ A s and a 
permutation matrix P CT such that 



We measure the size of the curvature vector v = (bi, 62 , &3, bi) of a Descartes configuration 



We claim that for positively oriented integral Descartes configurations with greatest common 
divisor g, we have 



and equality holds if and only if v is a permutation of (0, 0, g, g). If all curvatures are nonneg- 
ative this is clear, since at most two can be zero, and the other two are positive integers. Now 
in any Descartes configuration at most one circle can have negative curvature, call it 61 = —a 
(a £ in which case it encloses the other three. Each of these three enclosed circles has a 
larger curvature in absolute value than the bounding circle, so bi > a + 1 for i = 2, 3, 4. Thus 
size(v) > —a + 3(a + 1) > 2a + 3 > 2g, which proves the claim. 

We give a reduction procedure which chooses matrices in A s to reduce the size and show 
that the procedure halts only at a vector of form (0,0, g, g), up to a permutation. To specify 
it, we observe that for the curvature vector v = (b±, 62, 63, b<±) T of any integral Descartes 
configuration with &i < 02 < 03 < 04, we have 



P (7 U(6i,6 2 ,6 3 ^4) T = (0,0, 5 ,<7) T . 



by 




size(v) > 2g, 
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and equality holds if and only if bib 2 +b 2 b 3 = 0. To see this, we have S4V = (b\, b 2 ,b 3 , b' 4 ) T 
where 

64 = 2(6i + b 2 + 63) - h = h + 62 + h - 2y/bib 2 + b 2 b 3 + b 3 b 1 . 

Thus 

l T S 4 v - l T v = b' i -b A = -4 + Ms + 6361 < 0. 

Equality can hold if and only if b\b 2 + ^2^3 + ^3^1 = 0, which proves the observation. Note that 
b\ < when the equality in (|4.1j) holds. Also note that g = gcd(&i, b 2 , b$, 64) is an invariant 
under the action of A . 

Starting with any positively oriented integral Descartes configuration with curvature vector 
(61, 62, 63, 64) 71 , where 6, is the largest number, we apply S, G A By l|4.1jl . the size(v) decreases 
but cannot be negative, so after a finite series of Sj we arrive at a positively oriented integral 
Descartes configuration with curvatures v' = (b[, b' 2 , 6' 3 , b' 4 ) T , where gcd(6' 1 , b' 2 , 63, & 4 ) = 5 and 
the smallest curvature, say b[, satisfies b[ < 0, and the size of v' can not be reduced by the 
action of the Apollonian group A. Call this the basic reduction step. Note that the basic 
reduction step involves only matrices in the Apollonian group and therefore moves around 
inside a single Apollonian packing. 

If b[ = then necessarily b' 2 = 0, whence the curvature vector is (0, 0, b' 3 , b' 3 ) T , and by 
g = gcd(6' 1 , b 2 , b' 3 , b' A ), we have 6' 3 = g and the reduction halts. If b[ < 0, applying Sj , we get a 
new Descartes configuration with v" = (— 6^, b 2 + 2b[, b 3 + 2b[, b' A + 26^)^, which is positively 
oriented and lies in a new Apollonian packing and has 

size(y") = size(y) + 46' x < size(v). 

Thus the size strictly decreases and is non-negative. Now we may re-apply the basic reduction 
step. Continuing in this way we get strict decrease of size(v) at each step, with the only 
possible halting step being the smallest curvature equals 0. Since the size of the curvature 
vector is bounded below and decreases by at least one at each step, the procedure terminates 
at (0,0, g,g), up to a permutation. ■ 

Proof of Theorem 14. 1L For existence, the super-packing generated by a Descartes con- 
figuration with curvature vector (0,0, g, g), which is a homothetically scaled version of the 
configuration (0, 0, 1, 1) pictured in Figure 1, is necessarily integral with divisor g. 

For uniqueness, Theorem 14.21 shows that any two geometric integral Apollonian super- 
packings with divisor g each contain a Descartes configuration whose curvatures are (0, 0, g, g) 
up to permutation and orientation. Now it is true for any two such Descartes configurations 
with identical curvature vectors are congruent, i.e. one is obtainable from the other by a 
Euclidean motion. This is obvious by inspection for the (0, 0, g, g) Descartes configuration, 
which necessarily consists of two touching circles of radius | and two parallel lines. 

Now the Euclidean motion that takes one Descartes configuration to the other, also takes 
the super-packing generated by the first configuration to the one generated by the other, 
because the super-packing is defined by the action of the super-Apollonian group on the left 
on the Descartes configuration W©, and this commutes with the Euclidean motion acting as 
a Mobius transformation on the right. This establishes uniqueness. ■ 

We can use the freedom of a Euclidean motion allowed in Theorem 14. II to make an internal 
Apollonian super-packing strongly integral. 



13 



Theorem 4.3. For each integral geometric Apollonian super-packing there is a Euclidean 
motion that takes it to a strongly integral geometric Apollonian super-packing. 



Proof. Using Theorem 14.21 each integral geometric Apollonian packing contains a Descartes 
configuration D with curvatures (0,0, g,g); note that for a geometric packing the order and 
orientation of the Descartes configuration do not matter. The curvature vector determines 
the Descartes configuration up to congruence. We can now find a strongly integral Descartes 
configuration T>' with this curvature vector. For g = 1 such a configuration is given explicitly 
by 1)6. lj) below, and for larger g we obtain a strongly integral configuration from it using the 
homothety (x,y) i— ► -(x,y). There exists a Euclidean motion that maps T> to T>', since they 
are congruent configurations. This motion maps the super-packing to the super-packing 

,4 5 [2?'], which is strongly integral. ■ 



5. Strongly Integral Super-packings 

A Descartes configuration T> is strongly integral if its associated 4x3 curvature center- 
coordinate matrix is an integer matrix; this property is independent of ordering or orien- 
tation of the Descartes configuration. A super-packing is called strongly integral if it contains 
one (and hence all) Descartes configurations having this property. Since a strongly integral 
super-packing is integral, it has a divisor g as an invariant. 

Our main object in this section is to classify strongly integral super-packings, as follows. 

Theorem 5.1. (1) For each g > 1 there are exactly 8 different geometric Apollonian super- 
packings that are strongly integral and have divisor g. 

(2) The set of all ordered, oriented Descartes configurations that are strongly integral and 
have a given divisor g fill exactly 384 orbits of the super-Apollonian group. 

This theorem classifies these super-packings as rigid objects, not allowed to be moved by 
Euclidean motions. To prove this result we derive a normal form for a "super-root quadruple" 
in a super-packing of the kind above, as follows. 

Theorem 5.2. Given a strongly integral Apollonian super-packing A s [T>o] with the divisor 
g > 1, there exists a unique "reduced" Descartes configuration V G .A^fDo] whose curvature- 
center coordinate matrix M = is of the form A mtn [g] or B mjn [g] for m,n £ {0, 1}, up to 
a permutation of rows, where 









1 








-1 


9 


m 


n 


9 


m — 2 


n 



B OT ,ri [g] 



1 
0-1 

g m n 

g m n — 2 



(5.1) 



and the sign is determined by the orientation ofT>Q. 

Proof. For a strongly integral Apollonian super-packing .4 s [2?o] with the divisor g, by The- 
orem there exists a strongly integral Descartes configuration P £ .AsfPo] with curvatures 
i(0, 0, g,g). The two straight lines in T> must be parallel to either x-axis or y-axis. It follows 
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that the 4x3 curvature-center coordinate matrix is of the form A„ hn [g] , or B mjn [g] , for 
some m,fi£Z, up to a permutation of rows. 

We now reduce m, n to take the values 0, 1 using the following identities, which are easy to 
check. 

P(34) S 3 A m = A m _ 2) n[g], P( 34 )S3B min [g] = B mjn _ 2 [<7], 

P( 3 4)S4A min [<7] = A m+ 2 >n [g], P(34)S4B m>n [5(] = B min+ 2b]) 

F > (i2)Si A mt n[g\ = A m ^ n+ 2 [g] , P^Sj B m>n [^] = B m+ 2 )n [g], 

P(i2)S 2 r A miTt [5'] = A rn , n ^2[g], P(i2)S 2 r B m)n [c/] = B m _2, n [5']. 

where P(jj) is the permutation matrix that exchanges i and j. Also note that P<jSj = Sg-^Po-, 
Po-S^ = S^^P CT . Hence there is a series of group operations in A s which takes to 
a permutation of A mtn [g] or B mi „[<7] with m, n G {0,1}. This proves the existence of the 
"reduced" Descartes configuration in the Apollonian super-packing A 5 [Do]- 

To prove the uniqueness, it suffices to show that the 24x8x2 = 384 Descartes configurations 
whose curvature-center coordinate matrices are 

{PA m ,n[s],PB m ,n[j] I P G Perm 4 ,m,n G {0,1}} 

are in different Apollonian super-packings. (There are two signs for each of A mjn [g\ and 
~Bm,n[g]- ) In what follows we let A m ^ n [g], and B m)n [g] denote the unique 4x4 augmented 
curvature-center coordinate matrices extending A mi „[<?] and B m)n [p], respectively; uniqueness 
holds by Theorem 3.1 of part I. 

Note that each Sj and Sf preserves the (total) orientation of the Descartes configu- 
ration, as well as the parity of every element of Mp. First we show that the matrices 
A m>n [g] , B mi „ [g] , (m, n G {0,1} are in distinct orbits of A s x Perni4. To see this, for each 
integral vector veZ 4 , let k(v) be the number of even terms in v, and for any strongly integral 
Descartes configuration V let 

k(M v ) = (k(vi),k(v 2 ),k(v 3 )), 

where Vi,V2,V3 are the column vectors of Mp. Then n(Mp) is invariant under the action 
of A s x Perm4. For m,n G {0, 1}, K(A m:Tl [g]), K(B mjJl [gi]) are all distinct except k(Ai 5 o[<?]) = 
k(Bo,i[p]) = (*,2,2), where * is 2 if g is odd, and 4 if g is even. However, Aio[s] an d 
Bcitfi 1 ] can not be equivalent under the action of A s x Perm4. Arguing by contradiction, 
assume that there exists a matrix U G A s x Perni4 such that UA^ofg] = Bo,i[<?]. This 
relation lifts to augmented curvature-center coordinates: UA^ofg] = Bo 5 i[g]. It follows that 
U = (Ai i o[<7]) _1 Bo ) i[5'] is unique. We can directly verify that for m, n G Z 



^-m,n[g] 



&m,n[g] 



2(n + l)/g 
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2(1 -n)/g 
(m 2 + n 2 — l)/g 
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2(m + l)/g 





1 





2(1 -m)/g 
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(m 2 + n 2 - l)/g 
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n-2 



(5.2) 
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by checking that these satisfy the identity of Theorem 3.2 of part I necessary and sufficient to 
be augmented curvature-center coordinates. Now it is easy to verify that P(i4)P(23)DAi i o[g] = 

Bo,i[g], where D = — Qd is defined as in §3 of Part I [3], 



D 



•1 1 

1 -1 

1 1 

1 1 



1 1 

1 1 

-1 1 

1 -1 



(5.3) 



By the uniqueness P(i4)P(23)D 



U G A s x Pernio which is impossible since A s x Perni4 
consists of integral matrices only, while D has half integers. In conclusion, A m>n [g], B min [g], 
(m, n G {0, 1}) are in distinct orbits of A s x Perni4. 

The final step is to show that for any permutation P / I, PA mtn [g] (resp. PB mjn [g]) can 
not be obtained from A mj1l [g] (resp. B mjTl [g]) by an action of A s . That is, we claim: if for a 
permutation matrix P G Perni4, there exists a matrix U G A s such that 

UA m;n [g] = PAm in \g], or UB m>n \g] = PB miTl \g], 

then P = I. 

To establish the claim, consider again the 4x4 ACC-coordinate matrices A m ^ n \g\ and 
B mn [g], From §3.1 of Part I [3], for any Descartes configuration T>, the curvature-center 
coordinate matrix Mp can be uniquely extended to a 4 x 4 ACC-coordinate matrix Wp. It 
follows that if UA mj „[g] = PA mi „[j], then the equality holds for their 4x4 ACC-coordinate 
matrices, i.e., UA m ^[g] = PA mi „[g]. It implies U = P G A s n Perm4. However, comparing 
the size of U and P, where the size of a matrix U is defined as /(U) := 1 T U1, we have 
/(P) := 1 T P1 = 4 for P G Perm 4 , and /(U) := 1 T U1 > 8 for any U G A s , U/I, (c.f.§5 
of Part I Therefore the only possibility is U = P = I. The same argument applies to 
B m ,n[g], and the claim follows. 

We conclude that a reduced Descartes configuration of the form (|5.2j) in any strongly 
integral Apollonian super-packing exists and is unique. ■ 



Proof of Theorem 15. 1L (1) Since there are 48 orbits of ordered, oriented Descartes config- 
urations corresponding to each geometric super-packing, to show there are exactly 8 geometric 
super-packings it suffices to show that the strongly integral Descartes configurations form 384 
orbits of the Apollonian group, which we do below. 

(2) We enumerate the complete set of ordered, oriented Descartes configurations that are 
strongly integral, with greatest common divisor g, as follows. By Theorem 15.21 any such 
Descartes configuration is equivalent under the action of A s to a permutation of a Descartes 
configuration whose 4x3 curvature-center coordinate matrix is of the form A mtn or B mn , with 
m,n G {0,1}. The uniqueness of Theorem 15.21 asserts that the 24 permutations of ^4 mn [^] 
(B mtn [g]) are all in distinct orbits of the super- Apollonian group. Considering the two choices 
of orientations, we get 24 x 8 x 2 = 384 orbits. ■ 
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6. Primitive Strongly Integral Super-packings 



The strongly integral superpackings having a given curvature g.c.d. g are each obtainable from 
one with g = 1 by homothety. A homothety (x,y) \—* r(x,y) changes all curvatures by ~ while 
leaving (curvature) x (center) unchanged. Applying the homothety with r = ~ takes a strongly 
integral super-packing with g.c.d. g to one with g.c.d. equal to 1. 

We now consider strongly integral super-packings having g = 1, which we call primitive 
super-packings. Results for them carry over easily to those with divisor g > 1 by applying a 
homothety. Theorem 15.11 showed there are exactly 8 such packings. 

For convenience we will single out a particular one of them and term it the standard strongly 
integral super-packing. We choose this to be the super-packing generated by the ordered, 
oriented Descartes configuration having 



1 

-1 

1 
-1 



so that 



2 

2 

1 

1 



1 

-1 

1 
-1 



(6.1) 



This corresponds to a (0,0, 1, 1) Descartes quadruple, with the centers of the two circles lying 
along the x-axis and the circles touching at the origin (0, 0). The associated geometric integral 
super-packing is the one pictured in §3. Results proved below for the standard super-packing 
apply generally to all eight primitive integral super-packings, using the Euclidean motions 
mapping between them described after the proof of Theorem 16. II below. 

We first show that the geometric standard strongly-integral super-packing has a large group 
of symmetries, which form a crystallographic group of the plane. 

Theorem 6.1. The geometric standard strongly integral super-packing is invariant under the 
following Euclidean motions: 

(1) The lattice of translations (x, y) ^ (x + 2, y), and (x, y) i— > (x, y + 2), 

(2) The reflections (x,y) t— > (—x,y), and (x,y) i— ► (x,—y). 

The crystallographic group generated by these motions is the complete set of Euclidean motions 
leaving the geometric standard strongly integral super-packing invariant. 



Proof. The key fact used is that the action of the super- Apollonian group on Descartes con- 
figurations commutes with the action of Euclidean motions acting on Descartes configurations 
as Mobius transformations. This was shown in part I [5J Theorem 3.3(4)]. 

(1) There is a Descartes configuration corresponding to T>q shifted by 2 in the x-direction 
and the y-direction; call them T>q and T>q, respectively. These are given by the actions of S4 
and Sj~, respectively. Treating the x-shift first, we then have 

A s [D ]=A s [V* ]=A s [t x (v' )], 

in which i x : z z + 2 is the Euclidean motion translation by v = (2, 0) as in Appendix A of 
part I, and the ordered, oriented Descartes configuration V is a permutation of Vq. Then the 
geometric super-packing associated to ^[Pq] is therefore identical with that of .A^fDo]) and 
that of ^4 s [ta;(I?g)] translates it by (2,0). Thus the geometric packing is invariant under this 
translation. The argument for translation by (0, 2) is similar. 
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(2) This geometric Descartes configuration T>q is invariant under the reflections (x,y) \— * 
(—x,y), and (x,y) i— > (x,—y) viewed as Mobius transformations. The effect of these transfor- 
mations on the ordered, oriented Descartes configuration is to permute its rows. It follows as 
in (1) that the associated geometric Apollonian super-packings are identical. 

To see that these motions generate the full group of Euclidean motions leaving the super- 
packing invariant, we observe that the full group acts discontinuously on the plane. This is 
because the image of the (0, 0, 1, 1) configuration is either left fixed, or else it moves a distance 
of at least two in some direction, so that its circles do not overlap. Thus it must be contained 
in a crystallographic group whose translation subgroup is given by (1) above. Now the only 
possibilities are to extend the group by a subgroup of the finite point group of motions leaving 
(0,0) fixed (of order 8) leaving the lattice Z[(0, 2), (2, 0)] of translations invariant. Here (2) 
gives an extension of order 4. No larger extension occurs by observing that otherwise the image 
of the (0, 0, 1, 1) and (—1, 2, 2, 3) configurations at the origin would cross themselves. ■ 

One can now check that the 8 primitive geometric strongly integral super-packings given 
by Theorem 15 . 1 1 are obtained from the standard strongly integral super-packing by 8 cosets of 
the Euclidean motions (x,y) i— ► (x + 1,2/), (x,y) i— ► (x, y + 1) and (x, y) i— ► (y,x) with respect 
to the symmetries in Theorem 16.11 They are specified by the location and orientation of the 
(0, 0, 1, 1) configuration. 

Our next result shows that every primitive integral Descartes configuration with no curva- 
ture zero occurs inside the geometric standard strongly integral super-packing in a specified 
location. 

Theorem 6.2. In the geometric standard strongly integral super-packing, for each (unordered) 
primitive integral Descartes quadruple (a,b,c,d) except for (0,0,1,1), there exists a Descartes 
configuration having these curvatures, such that the center of its largest circle lies in the closed 
unit square {(x,y) : < x < 1, < y < 1}. The location of the center of this largest circle is 
unique. 

Proof. To establish existence, we first show that a Descartes configuration of the curvatures 
occurs somewhere inside the standard integral super-packing. This holds because the super- 
packing generated by such a configuration is an integral super-packing with divisor 1, which 
by Theorem 14.11 is unique up to a Euclidean motion. Thus the standard strongly integral 
super-packing must contain an isometric copy of it. Now that we have such a configuration 
inside the packing, we can use the translation symmetries in Theorem 16.11 to move it so that 
its largest circle has center inside the half-open square {(x,y) : —1 < x < 1,-1 < y < 1}. If 
we have —1 < x < then we apply the symmetry (x,y) t— > (—x,y), while if —1 < y < we 
apply (x,y) t— > (x, —y), as necessary. 

To establish uniqueness, we argue by contradiction. If uniqueness failed, there would exist 
a Euclidean motion taking one of these Descartes configurations to the other. Since a single 
Descartes configuration generates the entire super-packing, we conclude that the super-packing 
is left invariant under this extra Euclidean motion. Theorem 16. II described all such automor- 
phisms and all of them, except the identity, map every point in the interior of the unit square 
square strictly outside the square. This contradicts the assumption that the center of the 
largest circle of the first configuration is mapped to that of the second, when at least one of 
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these points is strictly inside the unit square. In the remaining cases where both centers lie on 
the boundary, one shows they are must lie the same boundary edge and that the automorphism 
leaves this edge fixed, so they are identical. Note that this argument shows in passing that the 
largest circle is unique, once (0,0, 1, 1) is excluded. ■ 

We come now to the main result of this section, which asserts that the geometric standard 
super-packing contains a copy of every integral Apollonian circle packing in a canonical way. 
The circles in the geometric standard super-packing can be foliated into a union of geometric 
Apollonian packings. 

Theorem 6.3. (1) Each circle in the standard super-packing with center inside the half-open 
unit square {(x,y) : < x < 1, < y < 1} is the exterior boundary circle of a unique primitive 
integral Apollonian circle packing contained in the geometric standard integral super-packing. 

(2) Every primitive integral Apollonian circle packing, except for the packing (0,0,1,1), 
occurs exactly once in this list. 

Proof. (1) Let the circle C be given. Recall from the proof of Theorem 13.11 that there is 
a unique minimal length admissible sequence U m U m _i • • • Ui[2?o] of generators of the super- 
Apollonian group that yield a Descartes configuration T> containing the given circle C, say in 
its j-th position. Admissibility requires that U m = Sj or Sf~ for some i ^ j. Then multiplying 
by Sj~ also gives an admissible sequence, and the Descartes configuration 

V := Sf[V] = S+U ro U TO _i • • • Ui[2>o] 

consists of the circle C plus three new circles nested inside the interior of C. This Descartes 
configuration T>' generates an Apollonian packing having the circle C as outer boundary, con- 
tained in the standard strongly integral super-packing. It is unique, because if there were a 
second Apollonian packing inside the bounding circle it would contain circles crossing those in 
the first packing, contradicting Theorem 13.11 

(2) Recall from Sect. 3 and 4] that inside each integer Apollonian packing is a positively 
oriented Descartes configuration whose absolute values of curvatures (a, b, c, d) are minimal, 
which is called a root quadruple. Theorem 4.1 of [3] showed that aside from the root quadruple 
(0,0,1,1) every root quadruple is of the form a<0<b<c<d. Root quadruples are 
characterized by satisfying the extra condition 

a + b + c>d. (6.2) 

All the circles in the resulting Apollonian packing are contained inside a bounding circle of 
curvature N = \a\, i.e., radius j^. Theorem 16.21 shows that for each root quadruple, with all 
curvatures nonzero there is a matching Descartes configuration whose largest circle has center 
inside the unit square and this largest circle is unique. The Apollonian packing contained 
inside this circle by is the integral packing with the given root quadruple, and it is unique by 
the result of (1). (In some cases, like (—1, 2, 2, 3) the root configuration is not unique, but the 
root quadruple and the packing itself are always unique.) Thus every primitive integer Apol- 
lonian circle packing, except (0,0, 1, 1), occurs exactly once bounding circles having center in 
the closed unit square. ■ 
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The initial part of the standard super-packing to depth 200 inside the unit square is pictured 
in Figure |H1 below. 

mod 1 




Figure 8: A "deeper" initial part of a super-packing (square of sidelength 1) 

One can make further computer experiments plotting the circles having various curvatures 
restricted (mod 4) inside the unit square. The results for circles having curvatures 1 (mod 2), 
2 (mod 4) and (mod 4) and size at most 200 are pictured in the following three figures. 

The figures empirically indicate that the following extra reflection symmetries occur; the 
line of symmetry is indicated on the figure with a dotted line. 

(a) 1 (mod 2) circles are symmetrical under reflection in the line x = 1 — y. 

(b) 2 (mod 4) circles are symmetrical under reflection in the horizontal line y = \. 
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Figure 9: Circles of curvature 1 (mod 2) 
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Figure 10: Circles of curvature 2 (mod 4) 
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(c) (mod 4) circles are symmetrical under reflection in the vertical line x = |. 

Based on this experimental evidence, one of the authors (CLM) conjectured that these 
symmetries hold. They were subsequently proved by S. Northshield |lUj . 

We may also illustrate these symmetry properties at the level of circles of a fixed curvature. 
Recall from §3 that each such circle contains a unique Apollonian packing having it as outer 
circle. Such a circle can then be labelled by root quadruple in the sense of [3] of this integral 
Apollonian packing. The cases of curvature —6, —8, and —9 corresponding to the three cases 
above were pictured in §2. In our paper considering number-theoretic properties of integral 
Apollonian packings, it was shown in |S[ Theorem 4.2] that the number of distinct primitive 
integral Apollonian packings with a given curvature — n of the outer circle had an interpreta- 
tion as a class number /i =t (— 4n 2 ) of positive definite binary quadratic forms of discriminant 
A = — 4re 2 , under GL(2, Z)-equivalence. One can raise the question whether there is some 
interpretation of the extra symmetries (a)-(c) above in terms of the associated class group 
structure under SL(2,Z) equivalence. 

7. The Super- Apollonian Group has Finite Covolume 

In this section we show that the super-Apollonian group is of finite volume as a discrete 
subgroup of the real Lie group Aut(Qx>) — 0(3, 1). This follows from the fact that the integer 
Lorentz group 0(3, 1; Z) has finite covolume in 0(3, 1), and the following result. 

Theorem 7.1. (1) The super- Apollonian group A s is a normal subgroup of index 48 in the 
group Aut(Q£ ) , r L). The group Aut(Qu,7j) is generated by the super-Apollonian group and the 
finite group of order 48 generated by the 4x4 permutation matrices and ±1. 

(2) The super- Apollonian group A is a normal subgroup of index 96 in the group G = 
Jo 0(3, 1; Z)J ~ 1 ; where 

1111" 
1 1-1-1 
1-1 1 -1 ' 
1-1-1 1 

The group G is generated by Aut(Q£ ) ,'L) and the duality matrix D, and consists of matrices 
with integer and half-integer entries. 

The duality matrix D is given by (|5.3j) in §5. This result allows us to associate to the 
super-Apollonian group the normal subgroup A s := Jq 1 A S 3 of the integer Lorentz group 
0(3, 1;Z), of index 96. 

The proof of Theorem 17.11 will be derived in a series of four lemmas. Let T be the group 
generated by adjoining to A the elements of the finite group of order 48 given by Perni4 and =LJ, 
and let f =< T,D >. The lemmas will prove that T = Aut(Qo,^) and f = Jo 0(3, 1;Z)Jq~ . 
Then analysis of cosets of A in these groups permits determining the index of A in these groups 
and showing normality. 

To determine T and T we will show 

T < Aut(Q D ,Z) <t = G. 




24 



It is easy to show that T is a subgroup of V of index 2, and D G T but D ^ Aut(Qr>, Z), so we get 
Aut(Q D ,Z) = T =< .A 5 , Perm 4 , ±1 >. It is easy to check the inclusion V < Aut(Qo,Z) since 
the generators Sj, Sf , P CT and ±1 of T are all in Aut(QD,Z). The inclusion Aut(QD,Z) < G is 
proved in Proposition 17,21 The equation G = V is proved in the next three lemmas. In Lemma 
17.31 we prove that the integer Lorentz group is exactly the group Aut(Lz) of invertible linear 
transformations that leave the integral Lorentz cone L% invariant. Lemma 17.41 states that in 
order to show that a group Q of invertible linear transformations of equals Aut(L%), it is 
enough to show that (1) the action of G on is transitive, and (2) there exists a point v G Li 
such that the stabilizer S v = {U G Aut(Lz) \ XJv = v} is a subset of Q. Using Lemma ITU we 
check that (1) the action of Jorj on L% is transitive, and (2), JoLJ contains the stabilizer S v 
of the point v = (1, 1, 0, 0) G L%. This proves G = T. 

Lemma 7.2. It is true that 



Aut(Q D ,Z) <G = J O(3, 1; Z)J^ l . 
Proof. Let U G Aut(Qr),Z). We need to show that the entries of Jq 1 UJq = JqUJq are all 



integers. Since Jq = ^H T — T, where 





11 

10 1 

110 



we have 



JoUJq = -(£ U v) llT - ^11 T UT - ^TU11 T + TUT, 



(7.1) 



h3 



where TUT is an integer matrix. 

From U T Q D U = Q D and Q D = |(2I - 11 T ), we have 

U T (2I - 11 T )U = 21 - 11 T . 



(7.2) 



Denote by Vj the i-th column of U, and by size{w ) := l T v the size of a vector v. Equating the 
entries of (|7.2j) we get 2vj • Vj — size(yi)size(y j) = 25ij — 1. In particular, size(vj), i = 1, 2,3, 4 
are odd integers. It follows that the matrix ^11 T UT is integral. 

Note that Aut(Qr>,Z) is closed under transposition. This is because 



U G Aut(Q D ,Z) U T Q D U = Q D 



U T Q D UQ D U T 

jT 



U 1 



=^ XJQdU 1 = {U 1 Qor'U 1 = (Qd)~ l = Qd^U 1 e Aut(Q D ,Z). (7.3) 

Applying the same argument of the preceding paragraph to U T , we then prove that the matrix 
±TU11 T is integral. 

Again using 2vj • v,- — size(vi)size(vj) = 25j,- — 1, summing over i, j = 1, . . . , 4, we get 



(E u *j 



2(vi + v 2 + v 3 + v 4 ) • (vi + v 2 + v 3 + v 4 ) + 8. 



i,3 



25 



Note that vi + V2 + V3 + V4 = (size(ri), size(r2), size{v^), size(ri)) where r-j is the i-th row 
of the matrix U. The sum ^ i (size(rj)) 2 is a multiple of 4 since each size(ri) is odd. It follows 
that j Uij) 2 is a multiple of 8. Since is an integer, we conclude that Ylij ls 

a multiple of 4 and the matrix Uij)ll T is integral. This proves Lemma 17.21 ■ 

Clearly A , Perm4 and ±1 are subgroups of Aut(Q£>,Z). Let T be the group generated by 
A s , Perm 4 and ±1, and let f := (r,D). Then T is a subgroup of T of index 2. 

The Lorentz light cone is the set of points {(yo> 2/1)2/2, 2/3 ) T G ^ 4 : — 2/o + 2/i + 2/1 + 2/3 = 0}- 
Let be the set of integer points in the Lorentz light cone, i,e, 

Lz := {(2/0, 2/1, 2/2, 2/3 ) T G Z 4 : - yl + y\ + y\ + y\ = 0}, 

and let Aut{L-£) be the set of linear transformations that leave invariant. 

Lemma 7.3. Aut(L z ) = 0(3, 1, Z). 

Proof. It is clear that 0(3, 1; Z) C ^4ui(L^). To show the other direction, let U G Aut(L%). For 
any integer point v G Lg, ?7v G Lg. Therefore (Uv) T Q,c(Uv) = 0, i.e., v T (U T Q,cU)v = 0. 
It is easy to check that the only symmetric matrices Q satisfying v T Qv = for all v £ L% are 
of the form Q = diag[— a, a, a, a]. Hence U T Q £ U = c 2 Q c , where c = det(U). 
Let 

"1 1 11" 
1-10 
1 ' 
1 _ 

Note that every column vector of X is an integer point in L%. Let Y = UX, and Z = U _1 X. 
Since U G Aut(Lz), every column of Y and Z is also an integer point in Therefore 
det(Y) = cdet(X) = —2c, and cdet(Z) = det(X) = —2. For each (2/0, 2/1, 2/2, 2/3) T G L%, we 
have 2/0 = 2/i + 2/i + 2/| • It follows that either yo, 2/1, 2/2, 2/3 are an even, or yo and exactly one 
of yi, 2/2, 2/3 are even. In both cases, det(Y) and det(Z) are even. This forces c = ±1. Hence 
U T Q £ U = Q £ , i.e., UeO(3,l). 

Since U maps points (1, ±1, 0, 0) T , (1, 0, ±1, 0) T , (1, 0, 0, ±1) T to integer points, if (a, b, c, d) 
is a row of U, then a ± b, a db c, a ± d G Z. Thus there exist integers a', b' , c' , d! of the same 
parity such that a = %, b = |-, c = |-, d = y. However, by Equation ()7.3() UQ,cU T = Q^. 
This implies -a 2 + 6 2 + c 2 + d 2 = ±1. Therefore b' 2 + c' 2 + d' 2 = a' 2 (mod 4). Hence a', 6', c', 
must all be even, which means a, b, c, <i are integers, and U G 0(3, 1; Z). ■ 

Lemma 7.4. Let Q be a group of linear transformations that preserves L%. If the action is 
transitive, and there exists a point v G L% such that the stabilizer S v = {U G Aut(L%)\\Jv = 
v} C Q, then Q = Aut(Lz). 

Proof. Clearly Q C Aut(Lz). For any P G Aut(Lz), assume P(w) = v' . Since Q acts 
transitively on Lz, there exists G G Q such that G(v') = v. That is, GP(v) = v. So 
GP G S v C g, and then P G G -1 ^ = ^. This proves Aut(L z ) Cg. ■ 

Lemma 7.5. f = G = J O(3, 1; Z) Jq 1 . 
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Proof. By Lemma 17.31 it is sufficient to prove that 

J x f J = J f J = Aut(L z ). 



(7.4) 



It is straightforward to check that JoSjJo, JoS^Jo, JoP CT Jo and JoDJo are integer matrices. 
In particular, 



JqSiJc 



-1 

-1 
-1 



-1 
-1 


-1 



-1 
-1 
-1 





J0P34J0 = P34 and JoDJo = diag[l, —1, —1, —1]. Therefore JorJo GJ 0(3, 1; Z) = Aut(L%). 

For any integer point (y , y 1} y 2 , 2/3) G ^z, -2/o + v\ + 2/1 + vl = implies 
2/0 + 2/1 + 2/2 + 2/3 = (mod 2). Hence Jo(2/o, 2/i) 2/2, 2/3) T is integral. It follows that 

Jo(iz) = a 2, 03 1 04) G Z 4 : (ai, 02, 03, 04) are curvatures of a Descartes configuration }. 

Then Theorem 14. 21 implies that Jol^Jo acts transitively on the integral Lorentz light cone Lz, 
where —I exchanges the total orientation of a point in Lz- 

We use Lemma l7~H to prove equation Q7.4JI with Q = JoIVo- Let v := (1, 1,0,0) G and 
consider its stabilizer 

S v := {U G 0(3, 1; Z) : Uv = v, U T Q £ U = Q £ }. 
Assume U = (ui,j)fj = i- Solving the equations 

U(1,1,0,0) T = (1,1,0,0), U T Q £ U = Q £ , 
we obtain the following linear and quadratic relations between the entries of U: 



(7.5) 



U12 = 
U22 



1 - un, 
z 1 - U21, 



U13 
U32 



U23, 
-U31, 



u 14 = n 24 , 
u 42 = -U41, 



and 



2,2 
n 33 + u 43 



2,2 
^34 + U M 



1. 



■"13 = «31«33 + "41«43, 



«33«34 + «43«44 = 0, 
""14 = U31U34 - U41U44, 



1 



"21 



(«31 +«4l)- 



It follows that the matrix U can be expressed as 



U 



/ 1 + t —t gm + hn km + In \ 

t 1 — t gm + hn km + In 

m —m g k 

\ n —n h I j 



(7.6) 



where t = (m 2 + n 2 )/2, g 2 + h? = k 2 + I 2 = 1, and gk + hi = 0. Since g, h,k,l G Z, we must 
have (g,h),(k,l) G {(±1,0), (0,±1)}. 
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We can classify the matrices of the form (|7.6|) into four types, up to a possible multiplication 
by P34) , as follows. 



Type I : 



/ 1 + 1 —t m n \ 

t 1 — t m n 

m —m 1 

n —n 1 / 



Type II : 



/ 1 + 1 —t m —n \ 

t 1 — t m —n 

m —m 1 

n —n —1 ) 



Type III : 



/ 1 + 1 —t —m n \ 

t 1 — t —m n 

m —m —1 

y n —n 1 / 



Type IV : 



/ l + t 


-t 


—m 


— n 


\ 


t 


1-t 


—m 


—n 




m 


—m 


-1 







\ n 


—n 





-1 


/ 



where t = (m 2 + n 2 ) /2 and m,n,t € Z. 

We denote a matrix of type X with parameters m, n by U(m, n; X). The following equations 
can be easily checked. 



U(m,n;I)U(fc,l 
U(m,n;n)U(Jb,Z 
U(m,n;m)U(Jb,Z 
U(m,n;IV)U(A:,/ 



I) = U(m + A;,n + / 

I) = U(m + k,n — I 

I) = U(m — k,n + I 

I) = U(m — k,n — I 



I) : 

II) , 

III) , 

IV) . 



Also we have 



UCM;!)- 1 =U(-1,-1;I), U(l, -1; I)" 1 = U(-l, 1; I). 



Therefore the stabilizer 5„ is generated by P 34 , U(0,0;I), U(0,0;II) U(0,0;III), U(0,0;IV) 
together with A = U(l, 1; I), B = U(l, -1; I). 

Note that U(0,0;I) = U(0,0;II) 2 , U(0,0;III) = P 34 U(0, 0; II)P 34 , and U(0,0;IV) = 
U(0,0;II)U(0,0;III). Thus S v is generated by the matrices P34, A, B, and C = U(0,0;II). 

Now one can check that 

P34 = J0P34J0 £ JorJcb 
C = diag[l, 1,1,-1] = J (-IP 2 3Pi4D)J G JofJo, 
A = (J SiJ )diag[l, 1,-1,-1] = Jo(S 1 Pi 2 P 34 )Jo e JofJo, 

B = (JoP34«Jo)C(JoS2Jo)C G JofJo- 



Hence S v C J rJo. By Lemma 17.41 JnFJn 
JqO(3, 1; Z)Jq 1 . This finishes the proof. ■ 



Aut(Lz 



0(3, 1;Z), or equivalently, T 



Proof of Theorem 17.11 Lemma 17.21 and Lemma 17.51 show that F < Aut(Qf),Z) < G. 
Hence V = Aut(Qu,7j) since T is a subgroup of G of index 2 and the duality matrix D is not 
in Aut(QD,%)- In other words, Auti y Qr ) ^l2) is generated by the super-Apollonian group A s 
and the finite group of order 48 generated by the 4x4 permutation matrices and ±7. The 



28 



super-Apollonian group A is a normal subgroup of Aut(QD,Z), since Po-SiP^-i = S CT (j), and 
P a SfP a -i = S^ (f j. The index is 48 since A s n (Perm 4 x {±1}) = I, (c.f. §5 of Part I [3]). 

By Lemma 17.51 the group G is generated by Aut(Qjj,Wi) and D. Note that D 2 = I and 
DSjD = Sf. It follows that the super- Apollonian group is a normal subgroup of G with index 
96. ■ 



The second part of Theorem 17.11 can be rephrased as asserting that the conjugate group 
A s = JqAJq is a normal subgroup of index 96 in 0(3,1). Its generators Sj = JoSjJg 1 and 
Sj- = JqS^Jq 1 are given by 



-1 

-1 
-1 



-1 -1 

-1 -1 

-1 

-1 



and £>2 



1 1 

1 1 

-1 

-1 



and 



1 

1 

-1 



-1 
1 

1 



1 

-1 
1 




and S4 



The generators Sj- = JqS^Jq 1 are given by 



(s,) r , 



1 -1 

-1 1 

1 

1 



which follows using Sj- = Sj and Jq = Jq 



8. Super-Integral Super-Packings 

This section treats the strongest form of integrality for super-packings, which is that where 
one (and hence all) Descartes configurations in the super-packing have an integral augmented 
curvature-center coordinate matrix Wp. We say that a Descartes configuration with this prop- 
erty is super-integral, and the same for the induced super-Apollonian packing. The following 
result classifies such packings. 

Theorem 8.1. (I) These are exactly 14 different geometric super-packings that are super- 
integral. 

(2) The set of ordered, oriented Descartes configurations that are super-integral comprise 
672 orbits of the super- Apollonian group. 

These packings are classified here as rigid objects, not movable by Euclidean motions. 
To prove this result, it suffices to determine which strongly-integral configurations are super- 
integral. The next result classifies the possible types of super-integral Descartes configurations, 
according to the allowed value of their divisors. 
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Theorem 8.2. Suppose that an ordered, oriented Descartes configuration T> in R 2 has integral 
curvature- center coordinates M = Mp, and let g = gcd(ai, 02, 03, 04), where (01,02,03,04) is 
its first column of signed curvatures. Then T> has integral augmented curvature- center coordi- 
nates Wj> if and only if one of the following conditions hold. 

(i) g = 1, or 

(ii) g = 2, and each row o/M has the sum of its last two entries being 1 (mod 2). 
(Hi) g = 4, and the last columns rows o/M are congruent to 



' 1 


" 




' 


1 " 


1 










1 


1 





or 





1 


_ 1 










1 _ 



(mod 2). 



Proof. By Theorem 15.21 there exists a matrix U £ A s and a permutation matrix P such that 

PUMp = A mtn \g] or B m>n [g], 

for some m,n £ {0,1}, where A mj „[g] and B m ^ n [g] are given in Formula ()5.1j) . and their 
corresponding augmented curvature-center coordinate matrices are ^4 mjn [o], B m<n [g], given in 
Formula (|5.2|) . Since each generator of A s preserves the super-integrality, as well as the parity 
of every element of Wp, it follows that Wj is integral if and only if one of the following 
conditions holds: 

1. g = 1, or 

2. g = 2, and m? + n 2 — 1 = (mod 2), or 

3. <? = 4, and the reduced form is ^0,1 [4] or £?io[4], up to a permutation of rows. 

In Case 2 we need m + n = l(mod2); in Case 3 the condition is equivalent to the one stated 
in the theorem. ■ 

Proof of Theorem EH (1) Since each geometric super-packing corresponds to 48 distinct 
orbits of the super-Apollonian group on ordered, oriented Descartes configurations, to show 
there are exactly 14 geometric super-packings, it suffices to show there are exactly 672 orbits 
of the super- Apollonian group that are super-integral, which is (2). 

f2)Theorems 15.21 and 18.21 allow us to classify the set of ordered, oriented Descartes config- 
urations that are super-integral by the action of A s . From the criterion of Theorem 18.21 we 
have: 

1. For o = l, any strongly integral Descartes configuration is super-integral. 

2. For g = 2, half of those orbits are super-integral, namely, those whose reduced forms are 
Aq^x [2] , Ai t o [2] , ^0,1 [2] or B lt o[2], up to a permutation of rows. 

3. For o = 4, one fourth of those orbits are super-integral, namely, those whose reduced 
form are Ao,i[4] or ^^[4], up to a permutation of rows. 

4. For g 7^ 1,2,4, there are no super-integral Descartes configurations. 
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g 


# of orbits of A 6 ' 


Representative 


(1,1,1, 


1) 


96 


Ai,i[l], B hl [l] 


(2,1,1, 


1) 


96 


Ai j0 [1], S ,i[l] 


(1, 1,2, 


1) 


48 


4),o[i] 


(1,1,1, 


2) 


48 


So,o[l] 


(4,1,2, 


1) 


48 




(4,1,1, 


2) 


48 


-Bi,o[l] 


(1,2,1, 


1) 


96 


Ai j0 [2], B ,i[2] 


(2,2,2, 


1) 


48 


4),i[2] 


(2,2,1, 


2) 


48 


5i,o[2] 


(1,4,2, 


1) 


48 


Ao,i[4] 


(1,4,1, 


2) 


48 


5i,o[4] 



Table 1: Orbits of super-integral Descartes configurations classified by g. 



More details of this calculation are given in Tabled To explain the notation in Tabled for 
any 4x4 integral matrix W, let be the greatest common divisor of entries u>i », W2,i, f 3,i> 
w^i in the i-th. column. Then the action of A s preserves the vector g = (51,52,53,54)- (Fot 
W = Wp, 52 is the greatest common divisor of the curvatures.) In each row of the table we 
give the number of orbits of A s formed by the set of ordered, oriented Descartes configurations 
that are super-integral with the given g. We also list the representatives of those orbits. Note 
that each entry in the column labelled "Representative" stands for 48 orbits, obtained by 
taking two choices of (total) orientation, and 24 choices of permutation of rows. Also note the 
symmetry that the configurations with g are inverse of the ones with g' = Pi2g- We conclude 
from the table that the set of ordered, oriented Descartes configurations that are super-integral 
comprise 384(1 + | + |) = 672 orbits of the super-Apollonian group. ■ 

9. Concluding remarks 

This paper showed that the ensemble of all primitive, strongly integral Apollonian circle pack- 
ings can be simultaneously described in terms of an orbit of a larger discrete group, the 
super-Apollonian group, acting on the standard strongly-integral super-packing. Study of the 
locations of the individual integer packings inside the standard super-packing, presented in §6, 
leads to interesting questions, not all of which are resolved. The standard super-packing also 
played a role in analyzing the structure of the super- Apollonian group as a discrete subgroup, 
carried out in §7. 

The various illustrations show the usefulness of graphical representations, as a guide to both 
finding and illustrating results. This contribution is due in large part to the statistician co- 
authors (CLM and AW). Graphics were particularly useful in finding extra symmetries of these 
objects, such as those illustrated in Figures 9-11 and subsequently proved by S. Northshield 
|10| . However one must not forget the adage of H. M. Stark |12| p. 225]: "A picture is worth 
a thousand words, provided one uses another thousand words to justify the picture." Section 
3 of this paper provides such a justification for certain features of Figure 4. 
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There remain some open questions, particularly concerning the classification of all integer 
root quadruples classified by fixed curvature — N of the bounding circle. This quantity is known 
to be interpretable as a class number, as described in Theorem 4.2]. In §6 we observed some 
symmetries of these root quadruples inside the standard super-packing, see Figures HHH1 There 
is a new invariant that can be associated to such quadruples, which is their nesting depth as 
defined in §4 with respect to the generating quadruple Vq of the standard super-packing. It 
would be interesting to see whether this invariant might give some further insight into class 
numbers. 
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10. Appendix. Strong Integrality Criterion 



This appendix establishes that to show a Descartes configuration is strongly integral it suffices 
to show that three of its four circles are strongly integral. This affirmatively answers a question 
posed to us by K. Stephenson. 

Theorem 10.1. A Descartes configuration T> has integral curvature- center coordinates 

if and only if it contains three circles having integer curvatures and whose curvaturex centers, 

viewed as complex numbers, lie in Z[i]. 

Proof. The condition is clearly necessary, and the problem is to show it is sufficient. We write 
the circle centers as complex numbers z,- = Xj + iyj. So suppose T> contains three circles with 
curvatures b\, b 2 , b 3 £i7L and with curvaturex centers &1Z1, b 2 z 2 , ^3 Z 3 £ Z[z]. We must show 
that the fourth circle in the configuration has 64 6 Z and 64Z4 6 Z[i]. 

For later use, we note that Theorem 3.1 of Part I [H] has a nice interpretation using complex 
numbers z to represent circle centers. This was formulated in |Hj as the Complex Descartes 
theorem. It gives 

b\z\ + b\z\ + b\z\ + b\z\ = -(61Z1 + 6 2 z 2 + 6 3 z 3 + 6 4 z 4 ) 2 . (10.1) 

and 

6?zi + b\z 2 + &3Z3 + b\z A = ^(Mi + b 2 z 2 + b 3 z 3 + & 4 z 4 )0i + b 2 + b 3 + 64) • (10.2) 

We claim that 64 £ Z. This is proved in the following two cases. 
Case 1. 616263 + 0. 

We first suppose that z\ = 0. If both x 2 and X3 are zero, then — ^- = j~ + ^, which 
means 6162 + 6263 + 6361 = 0. Hence 64 = 61 + 62 + 63 G Z. Otherwise by permuting Z2 and 
Z3 if necessary we may assume that x 2 7^ 0. Then the following equations encode the distance 
between the circle centers, since the circles touch. 

61 6 2 

4 + yl = (l i+ l/, (10.4) 

(a* " x 2 f + (y 3 - 2/ 2 ) 2 = + ^-) 2 . (10.5) 

62 63 

We wish to solve these equations for j/3 in terms of 61, 62, 63, cci and £2. To this end we subtract 
the first two equations from the third and obtain 

2(*2* 3 + vm) = (r + h 2 + (r + r) 2 - (r + r )2 := a 
61 6 2 61 63 6 2 63 

Calling the right side of this equation R, we obtain 

xz = — {R- 2y 2 y 3 ) . 
2x 2 
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where the division is allowed since X2 7^ 0. Substituting this in the second equation yields a 
quadratic equation in y 3 , with X3 eliminated, namely (after multiplying by 4x1), 



<4 + y 2 )y 2 - AR y 2 y 3 + R 2 - 4x 2 (l + i) 2 = 0. 

01 03 



Since y 3 is rational, this equation has rational solutions, so the discriminant A must be the 
square of a rational number. After some calculation one obtains 



Since x 2 ,bi,b 2 are nonzero it follows that b\b 2 + b 2 b 3 + ^i&3 is a perfect square. Viewing 
the Descartes equation as a quadratic equation in 64 we obtain the formula 64 := 61 + 62 + 

63 ± 2^/b\b 2 + &2^3 + bibs, which shows that both roots are integers. These are the oriented 
curvatures of the two possible choices for the fourth circle in the Descartes configuration, so 

64 6 Z. 

Now assume that z\ = x± + iy\ is arbitrary. Define 



Then s 2 ,t 2 , s 3>^3 are rational numbers, and they also satisfy equations ()10,3j) . (|1U.4|) . ()10.5|) . 
(just replace X2, 2/2, x 3 , 2/3 in l)lfl.3|) . (|10.4j) . (|10.5|) by s 2 , t 2 , S3, t 3 , respectively.) By the preceding 
argument, we again have 16s!(6a& 2 + + 6163) / (hhh) 2 = q 2 for some rational number q. 
This implies that b\b 2 + b 2 b 3 + b\b 3 is a perfect square and consequently that 64 is an integer. 

Case 2. b\b 2 b 3 = 0. This is proved similarly, with an easier calculation. 

The claim follows, so 64 £ Z. 

We now proceed to show that 64Z4 S Z[i]. Now gives 



The equation (|l().2j) gives 

(h + b 2 + b 3 - 6 4 )6 4 z 4 = 2(6 2 zi + &|z 2 + 63Z3) - (h + b 2 + b 3 + 6 4 )(Mi + b 2 z 2 + 6 3 z 3 ). (10.7) 

We treat two mutually exhaustive cases. 
Case 1. b\ + b 2 + 63 7^ 64. 

Then Q1U.7[) gives 64Z4 = x^ + iy^ for some rational numbers X4, 2/4. However b\Z\, b 2 z 2 , b 3 z 3 
are integers by hypothesis, whence (jlU.6f) shows that 64,24 is an algebraic integer. Since X4, ?/ 4 
are rational, we conclude that 64Z4 must be an integer, i.e., in Z[i]. 

Case 2. b\ + b 2 + b 3 = 64. In this case, we have b\b 2 + b 2 b 3 + b±b 3 = 0. Now from 1)10. 7j) . 




(s 2 ,t 2 ) := (x 2 - x 1 ,y 2 -y 1 ) and (s 3 , t 3 ) := (x 3 - xi, y 3 - y t ). 



64Z4 = &1Z1 + 6 2 Z2 + b 3 z 3 ± 2v / ^izi6 2 z2 + b 2 z 2 b 3 z 3 + 61Z163Z3. 



(10.6) 



we have 



> 2 zi + 6^2 + &3Z3 = (&i + &2 + &3X&1Z1 + b 2 z 2 + b 3 z 3 ) 



which can be simplified to 



b 2 b 3 zi + b\b 3 z 2 + b\b 2 z 3 = 0. 
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Thus 



(&izi)(6 2 z 2 ) + (hz 2 )(b 3 z 3 ) + (6 3 z 3 )(&izi) = b 1 (b 2 z 2 + b 3 z 3 ) blb3Z2 ^ &1 ^ 2Z3 + (b 2 z 2 )(b 3 z 3 ) 

—o 2 o 3 

= -blzl-bl4+z 2 z 3 (b 2 b 3 -^(bl + bl)) 
= -b\(z 2 -z 3 f. 

Now b\ (z 2 — z 3 ) is a Gaussian rational number whose square is integral. Hence bi {z 2 — z 3 ) 
must be integral. It follows that 64Z4 = b\Z\ + b 2 z 2 + 63Z3 ± 2b\(z 2 — z 3 )i is in Z[i]. ■ 

Since the Apollonian group consists of integer matrices, all Descartes configurations in 
an Apollonian packing generated by a strongly integral Descartes configuration are strongly 
integral. This explains the integrality properties of curvatures and curvature x center pictured 
in the packing in §1, for example. The previous result now gives a weaker necessary and 
sufficient condition for an Apollonian packing to be strongly integral. 

Theorem 10.2. An Apollonian circle packing is strongly integral if and only if it contains 
three mutually tangent circles which have integer curvature- center coordinates. 

Proof. Suppose we are give three mutually tangent circles in the packing that are strongly 
integral. Any set of three mutually tangent circles in the packing is part of some Descartes 
configuration in this packing. This follows from the recursive construction of the packing, which 
has a finite number of circles at each iteration. If iteration j is the first iteration at which all 
three tangent circles are present, at that iteration they necessarily belong to a unique Descartes 
configuration. Theorem 110.11 now implies that this Descartes configuration is strongly integral. 
It then follows that the whole Apollonian packing is strongly integral. This proves the "if" 
direction, and the "only if" direction is immediate. ■ 
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